1 



OV Chiral Symmetry and the Low-Energy Spectrum of the QCD Dirac 
Operator 

J.J.M. Verbaarschot, 
Department of Physics and Astronomy, 
Q i SUNY at Stony Brook, Stony Brook, NY 11794 

<N 

The order parameter of the chiral phase transition is directly related to the infrared 
part of the spectrum of the QCD Dirac operator. This part of the spectrum follows from 
the low energy limit of QCD which is given by a partition function of weakly interacting 
Goldstone modes. We find that the slope of the Dirac spectrum is determined by the 
^ ■ pion decay constant whereas for A <C 1/ L 2 Aq CD the correlations of the Dirac eigenvalues 
are given by a random matrix theory with the global symmetries of the QCD partition 
function. A possible observation of these continuum results in lattice QCD with staggered 
fermions is discussed. 
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P_i' 1- Introduction 

<D . 

It has been well established that QCD is the correct theory of the strong interactions. 
The vast body of evidence from lattice QCD simulations at low and intermediate energies 
^ ■ is complemented by perturbative calculations which become reliable at high energies 
where the renormalized coupling constant is small due to asymptotic freedom. In this 
lecture we will focus on another limit where QCD simplifies. Because of spontaneous 
breaking of chiral symmetry, QCD at low energy reduces to a theory of weakly interacting 
Goldstone bosons. Although this theory cannot be derived from QCD by means of an ab- 
initio calculation, its Lagrangian is determined uniquely by chiral symmetry and Lorentz 
invariance. The validity of this low-energy theory is based on the presence of a mass-gap 
which is a highly nontrivial and nonperturbative feature of QCD. 

One of the questions we wish to address is to what extent the existence of an effective 
low-energy theory imposes consistency conditions on the original theory. This question 
was first posed in for the mass dependence of the QCD partition function. They argued 
that for small quark masses it can be both obtained from the effective partition function 
and from QCD. Since the mass dependence of the QCD partition function is given by the 
average of the fermion determinant this imposes consistency conditions on the average 
properties of the eigenvalues of the QCD Dirac operator. 

Some of the properties of the Dirac eigenvalues are robust against large deformations of 
the gauge field action well outside the scaling window. This type of spectral universality 
has been investigated primarily within the context of Random Matrix Theory HQ . What 
has been found is that correlations of eigenvalues on the scale of the average level spacing 
are universal, i.e. they are robust against deformations of the probability distribution of 
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the matrix elements. The low-energy effective action is also insensitive to large deforma- 
tions of the gauge field action. The reason is the existence of a mass gap. In the next 
section we will relate this property to spectral universality. 

2. Valence Quarks and the Low-Energy Limit of QCD 

Because the same mass occurs both in the quark propagator and in the fermion deter- 
minant the average Dirac spectrum cannot be obtained directly from the QCD partition 
function. In order to access the spectrum of the Euclidean Dirac operator D one has to 
introduce the valence quark mass dependence of the chiral condensate defined by 

The average is over the Euclidean QCD action which includes a fermion determinant that 
depends on the the sea-quark masses. The spectral density per unit space-time volume 
V = L A of the Dirac operator D is directly related to to Y>(m v ), 

PW/V = ^ (E(*A + 6) - E(iA - e)) . (2) 
The generating function for S(m„) is defined by [§,[F]|§ 



In addition to the usual quarks this partition function contains valence quarks and its 
bosonic superpartners. The chiral condensates for the different (super-) flavors are given 
by the same expression in terms of the eigenvalues of the Dirac operator. We thus have 
a maximum breaking of the axial flavor symmetry. The low-energy modes are then given 
by the Goldstone modes associated with the spontaneous breaking of this symmetry. 
Their quark content can be either two sea-quarks, one sea quark and one valence quark, 
or two valence quarks. The low-energy effective partition follows from the flavor super- 
symmetry and its breaking and Lorentz invariance as is the case for the usual chiral 
Lagrangian |]8],|7||1- One major difference is that in this case the Goldstone manifold is a 
super- manifold with both compact and non-compact degrees of freedom |T0|,|7|,|9| . 

3. Scales in the Dirac Spectrum 

For a nonzero value of the chiral condensate E we can identify three important scales in 
the Dirac spectrum. The first scale is the smallest nonzero eigenvalue of the Dirac operator 
given by A min = l/p(0) = n/J^V. The second scale is the valence quark mass for which 
the Compton wavelength of the associated Goldstone bosons is equal to the size of the 
box. Using the Gell-Mann-Oakes-Renner relation we obtain as Thouless energy [pH| ,p^,^,|6| 



m c = — 2 , (4) 

where F the pion decay constant. A third scale is given by a typical hadronic mass scale. 
The three scales are ordered as A m i n « m c < A. For valence quark masses m v <C m c 
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the kinetic term in the effective action can be neglected and the low-energy partition 
function can be reduced to a zero dimensional integral. However, any partition function 
with a mass gap and the same pattern of chiral symmetry breaking as in QCD can be 
reduced this way. The simplest such theory is chiral Random Matrix Theory (chRMT) 



14[ . In that case spontaneous breaking of chiral symmetry arises in the limit of infinite 
matrices. The advantage of working with chRMT is that is relatively simple to derive the 
distribution of the smallest eigenvalues. Of course, the results for Y,(m v ) obtained from 
the partially quenched partition function and from chRMT coincide [0,0. 

The kinetic term is also determined uniquely by chiral symmetry and Lorentz invariance 
which allows us to calculate the Dirac spectrum in the domain m v <C A. This results in 
the slope of the Dirac spectrum at A = which for Nf massless flavors is given by [fFHT3f 

p'(0) _ ( Nf - 2)(N f + (3) Sp 

p(0) 1Qtt/3 F 4 ' { } 

Here, (3 denotes the Dyson index of the Dirac operator. For QCD with fundamental 
fermions and three or more colors (with (3 = 2) this result was first derived in [IS|. The 



other two possibilities, (3 = 1 and (3 = 4, refer to QCD with fundamental fermions and 
two colors and QCD with adjoint fermions and two or more colors, respectively. 

The domain below m c has been investigated extensively by means of lattice QCD sim- 
ulations and agreement with the chRMT results has been found [|5|JT~6| 23 1 . A somewhat 



surprising result is that the lattice QCD data reproduce the analytical result for zero 
topological charge. This will be explained in the next section. 

4. Approach to the Continuum Limit for Staggered Fermions 

The low-energy limit of QCD and the small Dirac eigenvalues are described by the 
same partition function. In order to recover the continuum Ua{1) symmetry of the stag- 
gered Dirac operator (without the Ua(1) symmetry) its smallest eigenvalues thus have to 
approach their continuum limit as well. The staggered Dirac operator can be written as 

D KS = D c + a 2 A 2 D R , (6) 

where Dq coincides with the continuum Dirac operator at low energies, a is the lattice 
spacing and A is a typical hadronic mass scale. The condition that the Dirac spectrum 
of Dxs approaches that of Dq is (with || • || the norm of an operator) 

\\a 2 A 2 D R \\ AXa. (7) 

With £ ~ A 3 , the spacing of the eigenvalues near zero in lattice units is given by AAa ~ 
l/p(0) ~ l/AfA d_1 a a! ~ 1 . Since ||£)r|| ~ 0(1), the condition (^) can be rewritten as 



a d+l \ d+l 



A d+1 < — or LA = N 1/d aA < jV*~*k. (8) 
N v 1 



But we also require a sufficiently large lattice with LA ^> 1 resulting in 

Nwm > i. (9) 
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Our naive estimate for the total number of lattice points for staggered fermions to ap- 
proach the continuum limit is given by iV (3 d+1 ) d . In two dimensions we need lattices of 
the order of 27 2 to be reasonably close to the continuum limit. This number is consistent 
with simulations of the Schwinger model with staggered fermions |T7|]. In four dimen- 
sions the situation is much worse. According to the same estimate continuum physics is 
only seen on lattices as large as 343 4 which explains that todays staggered lattices show 



agreement with chRMT results in the sector of zero topological charge [16,18,23,24 



5. Conclusions 

We have argued that the the distribution of the smallest eigenvalues of the Dirac op- 
erator is a signature for the pattern of chiral symmetry breaking of the QCD partition 
function. Both the correlations of the smallest eigenvalues and the slope of the Dirac 
spectrum from have been obtained from a chiral Lagrangian. The intercept of the Dirac 
spectrum determines the chiral condensate whereas its slope fixes the pion decay constant. 
However, it takes very large staggered lattices to reliably perform such analysis. 
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